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Abstract 

We investigate the influence of stellar migration caused by minor mergers (mass ratio from 1:70 
to 1:8) on the radial distribution of chemical abundances in the disks of Milky Way-like galaxies 
during the last four Gyr. A GPU-based pure N-body tree-code model without hydrodynamics and 
star formation was used. We computed a large set of mergers with different initial satellite masses, 
positions, and orbital velocities. We find that there is no significant metallicity change at any radius of 
the primary galaxy in the case of accretion of a low-mass satellite of 10® M© (mass ratio 1:70) except 
for the special case of prograde satellite motion in the disk plane of the host galaxy. The accretion of 
a satellite of a mass > 3 x 10® Mq (mass ratio 1:23) results in an appreciable increase of the chemical 
abundances at galactocentric distances larger than ^ 10 kpc. The radial abundance gradient flattens 
in the range of galactocentric distances from 5 to 15 kpc in the case of a merger with a satellite with 
a mass > 3 x 10® Mq. There is no significant change in the abundance gradient slope in the outer 
disk (from ~ 15 kpc up to 25 kpc) in any merger while the scatter in metallicities at a given radius 
significantly increases for most of the satellite’s initial masses/positions compared to the case of an 
isolated galaxy. This argues against attributing the break (flattening) of the abundance gradient near 
the optical radius observed in the extended disks of Milky Way-like galaxies only to merger-induced 
stellar migration. 

Subject headings: galaxies: abundances - galaxies: evolution -- galaxies: kinematics and dynamics 


1. INTRODUGTION 

It is well known that the heavy element abundance 
in the disks of the Galaxy and other spiral galax¬ 
ies decreases with increasing galactocentric distance. 
This radial metallicity gradient was revealed about 
50 years ago through measureme nts of the metallic¬ 
ity in stars and H II r egions (e.g.. lMavodlI976l : iSearlel 
Il971t iPeimbertI 119791 : [Shields fc Searlel 119781) . The 
oxygen abundance gradient in nearby spiral galaxies 
traced by H II region s ranges from 0 dex kpc~^ to 
^ —0.1 dex kpc~^ (e.g.. IVila-Costas fc Edmunddll992al 
Zaritskv et al.lHO^ iSanchez et al.l 12014 iPilvuein et al l 


2014f) . Estimations of the iron and oxygen abun¬ 


dance gradient in the Galaxy based on Cepheids 
and open clusters result in a slope of approx i¬ 
mately —0.05 dex kpc~^ (lAndrievskv et all l2002bllal: 


Bragaglia et al.ll^UUsI: iLemasle et al. 120081 : iPedicelli et alJ 


20091 : lYong et aHl2012h . However, the origin of the gra¬ 
dient and its relation to other macroscopic parameters of 
galaxies are still open for discussion. 

There are indications that radial metallicity profiles in 
the Milky Way and some giant spiral galaxies show a 
break within the optical radius i? 25 , which corresponds 
the isophotal radius at a B-band surface brightness of 
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25 mag arcsec“^. The iron abundance profile in the 
Galaxy traced by open clust ers is found to b e come shal- 
lower beyond 10~ 14 kpc by iBragaglia et al.l 1)200811 and 
lYong et Ml 1)201211 . The earlier data on Gepheids con¬ 
firm the flatt ening of the iron gradie nt in the outer 
Gala ctic disk (lAndrievskv et al.l l2002bl lal: ILemasle et al.l 
1200811 . while recent measurements of an extended sam¬ 
ple of Cepheids show less evident (if any) flattening of 
the gradient in the outer G alactic disk ()Luck &: Lamberd 
120111: ILemasle et aLll20I3ll . OB star observations also do 
not indicate an existenc e of a metallicity profil e break in 
the outer Galactic disk (iDaflon &: GunhalboOdH . The ra¬ 
dial distribution of gas phase oxygen and nitrogen abun¬ 
dances in the optical disks of the majority of spiral galax¬ 
ies traced by H II regions can be well fitted by a sin- 
gle exponential profi le (linear in the term of lg(X/H)) 
()Pilvugin et al.l[201^ while possible profil e slope changes 
are reported for some barred galaxies ([Martin &: RovI 
ll995HZahid fc Bresolinl[20Tll : iScarano et al.ll20Illl . 

Far-ultraviolet and near-ultraviolet observations by the 
Galaxy Evolution Explorer (GALEX) satellite reveal 
that a large fraction of late-type galaxies in the local 
universe show sig nificant star formation beyond their 
optic al radii R 25 (iGil de Paz et al.l 120071 : iThilker et al.l 
1200711 . Those galaxies are referred to as extended ultra¬ 
violet disk (XUV-disk) galaxies. Spectra of a number of 
H II regions in the extended disks of five spiral galaxies 
were measure d recently revealing flat gradients in the ex¬ 
tended disks (Bresolin et al.l 120091: IGoddard et al.l 1201 II: 


Bresolin et al. ^121: iPattemon et al.l 120121 : iWerk et al l 

2 OIIII . iPilvugin et al.l ( 201211 found the same radial gradi¬ 


ent breaks for nitrogen abundances in those galaxies and 
noted that the change in the gradient slope is more dis¬ 
tinct in the radial distribution of nitrogen than of oxygen 













































































2 


Zinchenko et al. 


abundances. Observations of OB stars in NGC 3621 also 
support the metallicity profile flattening in the XUV-disk 
obtai ned from the study of H II regions (|Kudritzki et al.l 
120141 . Thus, the abundance gradient is flatter in the ex¬ 
tended disks {R > R 25 ) as opposed to the optical disks 
( R < R 9 ,f,) in the case of a number of spiral galaxies. 

iBresolin et al.l (|2012h have advocated that it is unlikely 
that in situ star formation could have produced and dis¬ 
tributed enough heavy elements to enrich the interstellar 
medium in the extended d isks to the presently observed 
values. iWerk et all (j201fl[ l have considered three plau¬ 
sible mechanisms that may explain the relatively high 
metallicity of the extended disks: radial redistribution 
of centrally generated metals, strong galactic winds with 
subsequent fallback, and a past interaction leading to the 
accretion of enriched gas. The l ast mechanism, ho wever, 
cannot be the sole explanation (jWerk et al.ll2nilll . 

The transport of the heavy elements from the inner to 
the extended disks is also suggested as a possible explana¬ 
tion of flat rad i al abu ndance gradients (|Werk et al.ll2011l : 
IBresolin et al.l l20I2h . Particularly, it was shown that 
the presence of a central bar and spiral structure may 
lead to rather st rong radial migrat i on of stars and gas 
in ga la ctic disks (iFriedli et al.llI994: iSellwood fc Binneyl 
200^ iR oskar et al.l 1200^ iSchonrich fc Binneyl 200£; 
Minchev fc FamaevI I20I0I : iMinchev et al.l I20I3L I20I4 ) 
and a strong flattening during galaxy evolution in the 
[Fe/H] radial profile s for the old stellar populations 
(jMinchev et al.l [20141 . But while such a mechanism 
is quite effective in the inner disk, radial migration 
studies that are in line with or going back to the 
analytic prescription for r adial migration strengths of 
iSchbnrich fc Binii^ (|2009ll suggest very low migration 
for the outer disk in quiescent galaxies. 

Radial migration can be also caused by the interaction 
and merging with other galaxies. Cosmological simula¬ 
tions demonstrate that minor mergers should be a fre¬ 
quent and common process during the lifetime of galax¬ 
ies, in particular since z ^ I, and the number of satel¬ 
lites that fall on a Milky Way-like g alaxy strongly in¬ 
creas es with decreasing satellite mass (jKazantzidis et ahl 
1200811 . At the present day, we obserye the strong in- 
teractio n of the Milky Way w ith the Sagittarius dwarf 
galaxy (|Ibata et al.llI994 II995I1 . The center of this dwarf 
spheroidal (dSph) galaxy with a mass of < 10® Mq 
is located nearly behind our Gal actic center at a he - 
liocentric distance of ~ 26 kpc (|Monaco et al.l 120041 . 
A satellite with the mass of a few times 10® M 0 
can significantly heat the outer Galactic disk, excite 
the central spiral structure and create a warp that 
can i nduce significant streams in the yelocity distribu ¬ 
tion (iShumakoya fc BerczIS 120051 : iQuillen et al.1120091 1. 
iBird et ^ (|2012[ 1 considered the dynamical eyolution of 
a Milky Way-like galaxy with an accretion history mo- 
tiyated by cosmological simulations and concluded that 
satellite perturbations appear to be a distinct mecha¬ 
nism for inducing radial migration. Thus one may ex¬ 
pect a considerable role of minor mergers in the cheino- 
dynamical eyolution of galaxies during the last few Gyr. 

So far the attention has been mainly focused on the 
behayior of the radial metallicity gradient during the 
last few Gyr under the “major merg er” scenario, i.e., in- 
teractions between massiye galaxies (iKewley et al.ll201^ 


Montuori et al.l I20I0I : iRuoke et al.1 I20I0I : iPerez et’aTI 
201 If) , while less attention has been paid to the possible 
influence of stellar migration caused by a minor merger 
on the radial distribution of the metallicity in the disks 
of spiral galaxies. In this study we will carry out high 
resolution N-body/TREE-GPU simulations aimed at ex¬ 
amining the influence of pure stellar migration caused by 
minor mergers (with a mass ratio of the mergers from 
1:8 to 1:70) on the radial distribution of metallicity in 
the disks of large spiral galaxies out to 25 kpc. In par¬ 
ticular, we will inyestigate the role of the initial orbit 
orientation and of the mass of the merged satellite in the 
formation of present-day metallicity gradients. We will 
explore whether satellite accretion can plausibly be re¬ 
sponsible for the flattening of the metallicity gradient in 
Milky Way-like galaxy disks. 

2. MODEL AND CODE DESCRIPTION 

Our aim was to deyelop a simple model that can tell 
us how stellar migration caused by mergers may haye 
affected the radial distribution of metallicity in the disks 
of a Milky Way-like galaxy. In order to see clearly the 
purely dynamical mixing effect (migration) and to reduce 
the number of free parameters in our model we adopt 
a number of simplifications: (I) our merger system is 
purely based on N-body simulations, without gas; (2) 
we do not assume star formation; (3) the initial dynamic 
and kinematic parameters of a Milky Way-like galaxy are 
assumed to be close to the Milky Way parameters at the 
present time; (4) the initial metallicity distribution of our 
model is based on the measured abundance distributions 
in spiral galaxies traced by H II regions and young stellar 
populations; and (5) we neglect the possible dilution of 
disk metallicity by satellite stars deposited in the disk of 
the host galaxy. In this sense the satellite is considered 
just as a dark matter sub-halo. 

We consider galaxy eyolution only during the last four 
Gyr since the initial conditions of our model are close 
to those at the present time. Cosmological simulations 
of the formation of a Milky Way-like galaxy support our 
assumption that the disk of a Milky Way-like galaxy re¬ 
main ed quite stable during the last few Gyr (jBird et al.1 
I 2 OI 3 II . This approach leads us to a model with a strongly 
reduced number of degrees of freedom for the inyesti- 
gation of the influence of the migration of collisionless 
particles caused by minor mergers. 

For our Milky Way-li ke galaxy we adopt the 
iKiiiiken fc Dubinskil (Fl^ “MW-D” model of a disk 
galaxy with some slightly changed parameters. In our 
model we represent the galaxy with 2 million disk par¬ 
ticles with a total mass of 5.75 x 10^® Mq, 0.4 million 
bulge particles with a total mass of 1.29 x 10^® Mq and 
5 million of dark matter halo particles with a total mass 
of 9.89 X 10^^ Mq. We take the exponential scale length 
of the disk to be Rd = 2.97 kpc, the disk scale height Zd 
= 0.225 kpc, and the central radial yelocity dispersion ar 
= 123 km s“^. With this choice of parameters the stel¬ 
lar surface density at the solar yicinit y is ^^57 M^t^pc "® 
and Toomre’s stability parameter Q (jToomrel Il964f l is 
2.6. We report the total list of our host galaxy model 
parameters in Table [T] and Table [H 

The initial distribution of the particles in the satel¬ 
lite galaxy (“perturber”) is giyen by a simple spherically 
symmetric Plummer model with N^at = 33000 — 300000 
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Table 1 

Galaxy Model Parameters. 


bulge cut-ott potential, M/c 
Bulge velocity dispersion, crt 
Bulge central density, pb 

-4.7 

1.29 

125 

Halo central potential, 4^0 

-6.7 

Halo velocity dispersion, ao 

1.45 

Halo potential flattening, q 

1.0 

Halo concentration, C 

0.1 

Characteristic halo radius, i?a 

0.7 

Disk mass, 

1.1 

Disk scale radius, 

0.66 

Disk truncation radius, .Rout 

7.0 

Disk truncation width, SRout 

0.5 

Disk scale height, z^. 

0.05 


N ote. — All values of t he dim ensionless model were adopted 
from ll^iiken &: Dubinsl^ 1)199511 . The units of length, velocity, 
and mass used here are: RU = 4.5 kpc, VU = 220 km s~^, 

MU = 5.1 X 10^*^ M 0 , respectively. 


Table 2 

Galaxy Component Parameters. 


Number of bulge particles, TVb 
Bulge mass, Mb 

Bulge particle mass, mb 

4x10^ 
1.29x101“ Mq 
3.23x10"! M0 

Number of halo particles, 

5x10® 

Halo mass, Mb 

9.89x10“ Mq 

Halo particle mass, mb 

19.78x10'! Mq 

Number of disk particles, A^d 

2x10® 

Disk mass, Md 

5.75x101“ Mq 

Disk particle mass, md 

2.88x10'! Mq 

Disk scale radius, R^ 

2.97 kpc 

Disk scale height, 

225 pc 


particles and with an initial total mass from 10® M© 
{Nsat = 33000) to 9 X 10® Mq {Nsat = 300000). The 
initial half-mass radius was chosen to correspond to one 
third of the power of the initial mass of the satellite from 
.Rhm = 0.15 kpc to 0.31 kpc. In this way we keep the 
initial ratio of the tidal radius to the half-mass radius of 
the satellite galaxy roughly constant (in our case it was 
« 25). 

Oxygen was chosen as “metallicity indicator”. We 
adopted as an initial condition that all the particles at a 
given galactocentric distance have the same metallicity, 
which was taken from the radial abundance gradient re¬ 
lation. We assume the central oxygen abundance to be 
12-|-lg(0/H) = 8.8 and the slope of the gradient to be 
—0.05 dex kpc“^, i.e., 

12-f lg(§) =8.8-0.05-i? (1) 

H 


where R is galactocentric distance in kpc for the 
disk particles. These values are based on the mea- 
sured abundance distributions in spiral galaxie s (e.g. , 


Vila-Costas & EdmundsI I1992bl IZaritskv et alJ 

1994 

Maciel & Costal 201(1: IPilvugin et all 120031 12004 

2014 

Sanchez et al.ll2014D. This simplification seems to be ius- 


tified since detailed studies of the age-metallicity relation 
of stars in the Galactic disk (solar vicinity) show an al¬ 
most constant metallicity for the last four Gyr, accompa - 
nied by substantial scatter (e.g., iHavwood et al.l (I2013D . 


their figure 9: iCasagrande et ^ 1)20111) . their figure 18). 

Thus, on the one hand, our particles can be considered 
as “stars” since an N-body simulation is used and the 
metallicity of the particle does not change during the 
evolution. On the other hand, they can be considered as 
“gas clouds” since the metallicities of all the particles at a 
given radius at the initial time are taken to be the same. 
It should be noted that since our simulations deal with 
a pure N-body problem without involving any chemical 
enrichment, our results can be applied to other elements 
as well, not only to oxygen. There is only one restriction: 
the process of chemical radial mixing must be done by 
objects that can represented as a “N-body particles”. 

2.1. Code description 

We ran all simulations using our version of the 
ybody TREE- GRAPE (GRAvity PipEline) code 
l|Fukushige et al.l[20051 section 4.1). We also ported the 
original GRAPE-based tree code to many different hard¬ 
ware platforms including the GPU (with multi-thread 
usage under the SSE or AVX instructions) and also to 
the recent NVIDIA Graphics Processing Unit (GPU) 
platform using the Compute Unified Device Architecture 
(CUDA). In order to efficiently use the current acceler¬ 
ation hardware (GPU), we emplo yed a modified tree a l- 
gorithm, originally developed bv iBarnes fc Hud (1198611 . 
whi ch was first used on G RAPEAA bv iMakina (| 199 If ) 
and iFukushige et all 1)19911) . 

On a typical desktop hardware (CPU: i5-2500K with 
4 cores @ 3.3 GHz -|- GPU: GeForce GTX 570 with 
480 cores @ 1.46 GHz) we get the results for the 
full self-gravity force calculation routine with the typi- 
cal tree-construc t ion p arameters {9 = 0.5, Ug = 3500, 
()Fukushige et al.11200^ section 4.1.1^ for V = 1 M 
particles with the initial Plummer distribution in Ri 2 
sec. Our speed is quite comparable to the most ad¬ 
vanced an d recent fully GPU t ree code implementation 
(bonsai2: iBedorf et al.l ()2012allbl )). which does one full- 
force calculation on the same hardware for the same 
particle distribution and with the same opening angle 
in Ri 1 sec. In the last years we already used and ex¬ 
tensively tested our hardware-accelerator-based gravity 
calculation routine in a few of our galactic dynamics 
projects a nd get quite accurate results with a good per¬ 
formance (IShumakova fc BerczikI 120051: IBien et ^120081 : 
IPasetto et alJI20inil2011D' 

The current set of simulations was carried out with the 
GPU version of the code using local GPU clusters avail¬ 
able at the authors’ institutions (ARI: kepler, MAO: 
golowood, NAOG: laohu). 

We call the current version of our N-body code 
ber-galO. The source code of the program is publicly 
available from the FTP sit(0 of one of the authors. In 
our production runs for the gravitational force calcula¬ 
tion we use a very conservative opening angle 9 = 0.4. A 
gravitational softening parameter e = 40 pc was adopted 
for all the particles. We use the simple leap-frog integra¬ 
tion scheme with a fixed time-step At = 0.313 Myr to 
advance the particle positions and velocities during the 
calculation. One basic run covering up to 4 Gyr with 

7.4 million particles takes 55 hours. The total num¬ 
ber of models presented in the current paper is 34, so the 

® ftp://ftp.mao.kiev.ua/pub/users/berczik/ber-galO/ 
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total computational time was ^78 days. 

2.2. Disk stability 

First of all, we computed the evolution of the host 
(Milky Way-like) galaxy without a satellite (no “per- 
turber”) in order to test the stability of the structure and 
the radial oxygen abundance distribution during 4 Gyr of 
evolution. The first and second columns of Figure[T]show 
the column-density evolution of the disk-|-bulge compo¬ 
nents of the galaxy model at t = 0, 0.4, 0.8, 1.2, and 
4 Gyr. The face-on column density maps (left column of 
Figure[T|) show that the isolated galaxy forms a long-lived 
bar. 

Ghanges in the disk’s radial surface density, Toomre’s 
stability parameter Q, and in the oxygen abundance pro¬ 
files during the 4 Gyr of isolated host galaxy evolution are 
presented in the top panels of Figure [2j The radial sur¬ 
face density profiles of the disk show that the formation 
of a bar leads to radial mixing in the inner part of galac¬ 
tic disk. According to Toomre’s stability parameter Q, 
our disk model is stable against non-axisymmetric per¬ 
turbations (see Fourier analysis presented on Figure [TS]) 
with a minimal value Q = 1.9 at a radius of 3.6 kpc 
for the initial isolated disk setup with an increase of the 
minimal Q value during the run. 

Nevertheless, the variations of the disk’s radial surface 
density and radial metallicity profile during 4 Gyr of evo¬ 
lution of an undisturbed galaxy are small (top panels of 
Figure [2|). The maximum value of changes in the oxy¬ 
gen abundance across the entire range of galactocentric 
distances also does not exceed ^0.1 dex during 4 Gyr 
of evolution of our undisturbed galaxy. This can be con¬ 
sidered as evidence that the host galaxy model is suffi¬ 
ciently stable and can be used to study changes in the 
radial abundance distribution due to migration caused 
by the merging (interaction) with a satellite. 

2.3. Stability against disk and code parameters 

The top panels of Figure |3] show the influence of the 
variation of the disk parameters on the disk’s radial sur¬ 
face density (left column), Toomre’s stability parame¬ 
ter Q (middle column), and oxygen abundance profiles 
(right column). Profiles for model G02 (a galaxy with 
a basic disk scale radius Rd = 3 kpc and disk trunca¬ 
tion radius Rout = 31.5 kpc, which merged with a satel¬ 
lite of a mass of 6 x 10® Mq (black lines), and models 
with Rd = 3.9 kpc (green lines) and Rout = 25 kpc (red 
lines) are presented. The dotted lines indicate the ini¬ 
tial profiles, and the solid lines show profiles at 4 Gyr 
of evolution. The initial slopes of the oxygen abundance 
gradient are in all cases -0.05 dex kpc“^. It is evident 
that increasing the disk scale radius by 30% increases 
also the surface density of the outer disk, which leads 
to a decrease of Toomre’s stability parameter Q in the 
outer disk, but makes it more stable against changes of 
the oxygen abundance induced by mergers. This behav¬ 
ior can be explained as follows. The level of the increase 
of the abundance in the outer disk depends on the ratio 
of the number of high-metallicity particles that migrated 
from the inner region of the disk to the number of parti¬ 
cles in the outer disk. Since increasing the Rd increases 
the number of non-inigrated particles (surface density) 
in the outer disk, the increase the oxygen abundance in 


the outer disk is less for a disk with higher Rd- A re¬ 
duction of the disk size from 31.5 kpc to 25 kpc shows 
no appreciable changes in the oxygen abundance profile 
within the reduced disk radius, i.e., up to 25 kpc. 

The bottom panels of Figure [3] show the stability of 
the disk’s radial surface density (left column), Toomre’s 
stability parameter Q (middle column), and the oxy¬ 
gen abundance profiles (right column) for different time 
steps in the simulation. Profiles for model G02 with 
At = 0.625 Myr (blue line). At = 0.313 Myr (red line), 
and At = 0.156 Myr (green line) are presented. The dot¬ 
ted lines indicate initial profiles. The solid lines indicate 
profiles at t = 4 Gyr. 

3. SATELLITE INITIAL ORBITS AND VELOCITY 

We ran four sets of simulations with satellite masses 
ranging from 10® Mq to 9 x 10® Mq (mass ratios of 
1:70 to 1:8). Each set consists of 7 runs with different 
initial positions of the satellite galaxy. The initial dis¬ 
tance from the center of the host galaxy was taken to be 
40 kpc for all the runs. The initial three-dimensional po¬ 
sitions and velocity directions were chosen in such a way 
that the satellite orbit inclination is equal to the fixed 
angles * = 0°, 30°, 60°, 90°, 120°, 150°, and 180°. The 
inclinations i > 90° correspond to a retrograde orbital 
motion of the satellite. The initial absolute value of the 
satellite velocity for our runs is V = 100 km s“^ (see dis¬ 
cussion below). The parameters of the simulations are 
summarized in Table 01 

To check how the orbital elements affect the radial 
abundance distribution we ran a set of five simulations 
with an initial satellite velocity from V = 0 km s“^ to 
200 km s“^, with a fixed initial distance from the galactic 
center of 40 kpc, with a fixed initial orbital inclination 
i = 30°, and with a fixed mass of Mgat = 6 x 10® Mq. 
The parameters of this set of simulations are summarized 
in the Tabled The results are discussed below. 

Eigure 0 ] shows the orbital evolution of satellites for 
this set of models. The top panel shows the evolution 
of distance from the host galaxy to the satellite in the 
disk plane. The bottom panel shows the evolution of 
satellite height above the host galaxy disk plane. The 
times needed for settling the satellite in the host galaxy 
center vary from ~ 0.7 Gyr (initial satellite velocity 
equals V = 0 km s“^) to ~ 3 Gyr (satellite velocity 
V = 200 km s-i). 

Eigure [5] shows the orbital evolution of satellites in the 
models A02, B02, C02, and D02, which have the same 
satellite initial position Rxy = 35 kpc and Rz = 20 kpc 
(initial orbital inclination i = 30°) and masses from 10® 
Mg to 9 X 10® Mq. The top panel shows the evolution of 
distance from the host galaxy to the satellite in the disk 
plane. The bottom panel shows the evolution of satellite 
height above the host galaxy disk plane. 

The two right columns of Figure 0 ] show 
disk-bbulge-l-satellite column density for the model 
with a satellite mass of 6 x 10® Mq (mass ratio 
1 :12), initial satellite position Rxy = 35 kpc and 
Rz = 20 kpc (initial orbital inclination i = 30°) and 
velocity V = 100 km s“^ at the times t = 0, 0.4, 0.8, 
1.2, and 4 Gyr. For comparison, the two left columns of 
Figure [T] show the column density for the run without 
satellite. The bottom right panel of Figure [2] illustrates 


On the influence of minor mergers on the radial abnndance gradient in the disks of the Milky Way-like galaxies 5 


20 


-- 0 


-10 

-20 


20 


o 10 
a 


- 0 

>H 


-10 

-20 


20 


u 10 
a 


- 0 


-10 

-20 


20 


-- 0 
>H 


-10 

-20 


20 


o 10 
a 


- 0 


-10 

-20 


- 20-10 0 10 20 
X [kpc] 



- 20-10 0 10 20 
Y [kpc] 


20 


-- 0 


-10 

-20 


20 


o 10 
a 


- 0 

>H 


-10 

-20 


20 


u 10 
a 


- 0 


-10 

-20 


20 


-- 0 
>H 


-10 

-20 


20 


o 10 
a 


- 0 


-10 

-20 


- 20-10 0 10 20 
X [kpc] 


- 20-10 0 10 20 
Y [kpc] 





Figure 1. Comparison of the log-scale column density image for disk-|-bulge-Ssatellite components in the model AGO (1st and 2nd column) 
and model C02 (3rd and 4th column). Top to bottom: density at £ = 0, 0.4, 0.8, 1.2, and 4 Gyr. 


the evolution of the radial metallicity distribution in the 
host galaxy during the merger with the satellite. The 
solid black line shows the initial metallicity distribution. 
The dotted lines mark the metallicity distributions at 
different times up to 4 Gyr. The top right panel of 
the Figure [2] shows the evolution of the metallicity 
distribution in the galaxy withont satellite. 

The evolution of the metallicity distribution in the 
model C02 can be divided into two stages. During the 
first stage, for np to i ~ 1.2 Gyr, the satellite moves 
along a semi-periodic orbit around the host galaxy, per¬ 
turbing the disk of the host galaxy. These perturba¬ 
tions induce the radial migration of stars and, therefore, 
result in changes in the metallicity distribution across 
the disk, increasing the abundances in its outer region 
{R > 10 kpc). During the second stage, the satellite 
approaches the host galaxy and passes through denser 
regions. The rates of loss of both angular momentum 
and kinetic energ y increase due to the intense dynamical 
friction (see, e.g.. lJnst et al.l[20TTI ). This effect moves the 


satellite towards the disk and then into the center of the 
host galaxy: in model G02, the surface density defined by 
the disk particles at the end of the simulation (t = 4 Gyr) 
is more than 100 times larger than the surface density de¬ 
fined by the satellite particles for all disk regions except 
the innermost region (R < 5 kpc). The lines on the bot¬ 
tom right panel of Figure [2] show the radial abundance 
distributions in the host galaxy at the times t = 1.2, 3, 
and 4 Gyr. There are no significant changes in the oxygen 
abundances in the inner part of the disk (R < 6 kpc), but 
the abundances at large galactocentric distances increase 
by ^ 0.2 dex. The slope of the abundance gradient at 
galactocentric distances 5 < i? < 15 kpc becomes flatter. 

Figure [7| shows the angular momentum distribution 
evolution for models that have the same satellite initial 
position Rxy = 35 kpc and Rz = 20 kpc (initial orbital 
inclination i = 30°) and mass 6 x 10^ Mq), but different 
values of the initial velocity. Both axes are in kpc units 
dne to the normalization of the angular momentum by 
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Figure 2. Example of the evolution of the disk radial surface density (left column), Toomre’s stability parameter Q (middle column) and 
oxygen abundance profiles (right column) during the 4 Gyr considered in our simulations. Top: galaxy without satellite. Bottom: galaxy 
with a merging satellite with a mass of 6 x 10^ Mq (mass ratio 1:12). The satellite’s orbital inclination is 2 = 30° (model C02). 








Figure 3. Top: Influence of the variation of the parameters of the disk on the disk’s radial surface density (left column), Toomre’s stability 
parameter Q (middle column), and oxygen abundance profiles (right column). Profiles for model C02 (a galaxy with a basic disk scale 
radius = 3 kpc and disk truncation radius Rout = 31.5 kpc) (black lines), and models with R^ = 3.9 kpc (green lines) and -Rout = 25 kpc 
(red lines) are presented. The dotted lines indicate initial profiles. The solid lines show profiles at £ = 4 Gyr. Bottom: Stability of the disk 
radial surface density (left column), Toomre’s stability parameter Q (middle column), and oxygen abundance profiles (right column) for 
different time steps. Profiles for model C02 with At = 0.625 Myr (blue line), At = 0.313 Myr (red line), and At = 0.156 Myr (green line) 
are presented. The dotted lines indicate initial profiles; solid lines indicate profiles at t = 4 Gyr. 


the circular velocity in the disks. All rows except the last 
represent the successive changes in angular momentum 
within the time interval of At = 0.5 Gyr. The last row 
represents the total cumulative changes of the angular 
momentum from 0 to 4 Gyr. 

Figure [S] shows the radial abundance profile at f = 
4 Gyr. The solid black line is the model without satel¬ 
lite. The dotted lines represent models with satellites 
of different initial velocities. One can see in Figure |4] 
and Figure [7] that all the models undergo merging by 


t = 4 Gyr, so we can compare the influence of the satel¬ 
lite’s initial velocity on the metallicity gradient as the 
result of the merging process. The largest influence on 
the metallicity gradient of the outer disk can be observed 
for the satellite initial orbital velocity of V > 100 km s“^. 
A further increase of the satellite’s initial orbital velocity 
does not significantly change the behavior of the radial 
oxygen abundance prohle of the outer part of the disk 
(Figure [6|). Thus, we choose the satellite’s initial orbital 
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Table 3 

Model runs with different satellite initial velocities and the resulting disk oxygen abundance gradients. 


Model 

Mass [x 10^Mq] 

V [km s""] 

Alg(0/H),n 

Alg(0/H)inner 

Alg(0/H)outer 

AOO 

0 

- 

-0.0453±0.0010 

-0.0387±0.0017 

-0.0554±0.0007 

EOl 

6 

0 

-0.0449±0.0007 

-0.0416±0.0008 

-0.0544±0.0010 

E02 

6 

50 

-0.0404±0.0008 

-0.0364±0.0007 

-0.0531±0.0006 

EOS 

6 

100 

-0.0326±0.0006 

-0.0265±0.0007 

-0.0426±0.0007 

E04 

6 

150 

-0.0300±0.0008 

-0.0231±0.0008 

-0.0460±0.0014 

E05 

6 

200 

-0.0298±0.0007 

-0.0223±0.0004 

-0.0439±0.0018 


Note. — First column: Model series number. Column 2: Satellite mass in units of solar masses. Column 3: Initial velocity 
of the satellite in km s~^. Column 4: Oxygen abundance gradients from fitting radial abundance profiles across the disk radius 
range of 7? = 0 to 25 kpc in dex kpc“^. Column 5: Oxygen abundance gradients from fitting radial abundance profiles in the 
inner disk (77 = 5 - 15 kpc) in dex kpc“^. Column 6: Oxygen abundance gradients from fitting radial abundance profiles in the 
outer disk (77 = 15 - 25 kpc) in dex kpc~^. The first row of this table shows the results for a host galaxy without a satellite. 




Figure 4. Orbital evolution of satellites in models that have the 
same satellite initial positions 77xy = 35 kpc and Rz = 20 kpc 
(initial orbital inclination i = 30°) and mass (6 x 10® Mq) but 
different initial velocity. Top panel: evolution of distance from the 
host galaxy center to the satellite in the disk plane. Bottom panel: 
evolution of satellite height above the host galaxy disk plane. 



Figure 6. Radial oxygen abundance profile at 4 Gyr. The solid 
black line is a model without satellite. The dotted lines indicate 
satellites with different initial orbital velocities (models EOl — E05). 

velocity V = 100 km as the basic setup for further 
investigation. 



t [Gyr] 



Figure 5. Orbital evolution of satellites in our models A02, B02, 
C02, D02, which have the same satellite initial positions 77xy = 
35 kpc and Rz = 20 kpc (initial orbital inclination i = 30°) and 
masses from 10® Mq to 9 x 10® Mq. Top panel: evolution of the 
distance from the host galaxy center to the satellite in the disk 
plane. Bottom panel: evolution of satellite height above the host 
galaxy disk plane. 


4. RESULTS & DISCUSSION 
4.1. Radial metallicity distribution 

Now we discuss the influence of the particle migration 
caused by the merger of satellite and host galaxy with 
different mass ratios and orbit orientations on the abun¬ 
dance distribution across the disk. 

Figure [5] shows the radial abundance distributions at 
t = 4 Gyr. The masses of the satellites are 10® Mq 
(upper left panel), 3x 10® Mq (upper right panel), 6x 10® 
Mq (bottom left panel), and 9 x 10® Mq (bottom right 
panel). The solid black line corresponds to the evolution 
of the model without satellite, the dotted lines show the 
host galaxy accreting satellites with different initial orbit 
inclinations. The comparison between the model without 
and with satellites shows that there are no significant 
changes (^0.1 dex) in the radial abundance distribution 
after 4 Gyr of evolution in the inner part of the disk 
of the host galaxy for all the satellite masses considered 
(Figure [8]). In the case of the smallest satellite mass 
(10® Mq) (top left panel), there are no significant changes 
in the radial abundance distribution across the whole 
disk except for the case when the orbital inclination of 
the satellite i = 0°. Accretion of more massive satellites 
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Figure 7. Angu[ar momentum distribution evolution for models that have the same satellite initial position -Rxy = 35 kpc and Rz = 20 kpc 
(initial orbital inclination i = 30°) and mass 6 X 10® Mq but different initial orbital velocities (different columns). Both axes are in kpc 
units (see the text). Top to bottom: radial mixing in time intervals of At = 0.5 Gyr at (0 - 0.5) Gyr, (0.5 - 1.0) Gyr, (1.0 - 1.5) Gyr, (1.5 
- 2.0) Gyr, (2.0 - 2.5) Gyr, (2.5 - 3.0) Gyr, (3.0 - 3.5) Gyr, and (3.5 - 4.0) Gyr. Last row: total between (0 - 4) Gyr. The color coding 
indicates the percentage of particles below some fixed levels (see color bar on the bottom right). 
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causes a significant increase (up to > 0.2 — 0.5 dex) in 
the oxygen abundances at galactocentric distances R > 
10 kpc, as illustrated in the top right and bottom panels 
of Figure m 

The accretion of a satellite with a mass of Mg^t ^ 
3 X 10® Mq or less with an orbit in the disk plane in 
prograde direction can lead to a significant increase of 
the metallicity in the outer regions of the host galaxy. It 
should be noted that the abundance redistribution due to 
merging strongly depends not only on the satellite mass 
but also on its initial orbital inclination. The maximum 
redistribution is caused by satellites with prograde orbits 
(i = 0°) in the disk plane of the galaxy. The minimum 
redistribution is caused by satellites with retrograde or¬ 
bits with an initial inclination i = 90° to 150°. 

Figure |9] shows the mean deviation around the radial 
abundance trend crig(o/H) (in radial cylindrical bins) at 
t = 4 Gyr for models without (solid black line) and with 
satellite (dotted lines). Models with the four satellite 
masses 10® Mq (mass ratio 1:70) (top left panel), 3 x 
10® Mq (mass ratio 1:23) (top right panel), 6 x 10® Mq 
( mass ratio 1:12) (bottom left panel), and 9 x 10® M© 
(mass ratio 1:8) (bottom right panel) are presented. The 
scatter in the abundances at a given radius caused by 
radial migration of the particles in the model without 
satellite increases with galactocentric radius up to i? 

15 kpc and decreases beyond this radius. The scatter 
in metallicities at a given radius in the outer disk (> 
15 kpc) significantly increases for most satellite initial 
masses/positions compared to the isolated galaxy case. 

Figure [To] shows the abundance distribution histogram 
for particles at f = 4 Gyr for the models AOO (with¬ 
out satellite) and A02, B02, G02, D02 (with satellites of 
different masses) at galactocentric radii of i? = 10 kpc 
(top panel) and R = 20 kpc (bottom panel). The verti¬ 
cal black lines show the initial abundance distributions 
at these radii. It is clear that the accretion of a satel¬ 
lite with a mass > 3 x 10® Mq significantly shifts the 
abundance distribution at both galactocentric distances 
to higher values and the shift is larger at larger radii. A 
small fraction of particles with high abundances from the 
central region migrates to the outer disk region even in 
the case of an isolated galaxy. These particles produce an 
asymmetric particle abundance distribution with a high 
abundance tail (see the solid black line at the bottom 
panel in Figure Ho}. Figure |TT] shows again the abun¬ 
dance distribution histogram for particles at 4 Gyr for 
galactocentric radii of 10 kpc (top panel) and 20 kpc 
(bottom panel) predicted by the models with a fixed ini¬ 
tial satellite orbit inclination (i = 30°) and mass (Mgat 
= 6 X 10® Mq) but with different initial velocities from 
F = 0 to 200 km s'h 

4.2. Radial metallicity distribution breaks 

Our simulations provide radial metallicity distributions 
in the outer disk (up to 30 kpc) for a Milky Way-like disk 
galaxy. It gives us an opportunity to test whether the 
shallow abundance gradient observed in the outer disks 
of some spiral galaxies can be attributed to star particle 
migration caused by merging. The part of the disk with 
R between 5 and 15 kpc will be referred to as the inner 
(optical) disk while the part of the disk between 15 and 
25 kpc will be referred to as the outer (extended) disk. 


The innermost region i? < 5 kpc is not discussed here. 
The region beyond i? = 25 kpc is also excluded from our 
discussion since the number of particles there is not high 
enough for reliable estimations of metallicities. 

A close examination of the radial metallicity distri¬ 
butions shows that the abundance distribution slopes 
slightly change with radius. To quantify the variations 
in the abundance gradient slopes we fit the abundance 
distributions at 4 Gyr with separate exponentials in the 
inner and outer disks for each simulation from Table [3] 
and Table m The obtained slopes are listed in Table [H 
The first column gives the model number. Column 2 re¬ 
ports the satellite mass, column 3 lists the initial distance 
of the satellite, column 4 reports the initial satellite incli¬ 
nation, columns 5, 6, and 7 give the slope of abundance 
gradient in entire disk (i? = 0 to 25 kpc), in the inner 
disk (i? = 5 to 15 kpc), and in the outer disk (i? = 15 to 
25 kpc), respectively. 

Inspection of Table HI shows that the slopes of 
the abundance gradient at t = 4 Gyr in the in¬ 
ner disk (—0.039 dex kpc“^) and in the outer disk 
(—0.055 dex kpc“^) of an isolated galaxy (model AOO) 
are quite close to the initial slope of —0.05 dex kpc“^. 
The slopes of the gradient in the outer (extended) disk 
in the majority of the models with satellites are close to 
that in the model of an isolated galaxy. In contrast, the 
slopes of the gradient in the inner disk in the majority 
of the models with satellites of masses > 3 x 10® Mq 
decrease. Assumi ng a mass-to-light ratio in the B band 
of M/Lb = 1-4 (jFlvnn et al.ll200fil l. the optical {R 25 ) 
radius of our host galaxy is 15 kpc. Thus, the flatten¬ 
ing of the gradient observed in extende d ultraviolet disks 
(XUV-disk) of some spi ral galaxies (|Gil de Paz et al.l 
l2007t iThilker et al.l 120071) cannot be attributed to par¬ 
ticle migration caused by the accretion of a low-mass 
satellite at least for Milky Way-like type galaxies. 

4.3. Particle migration 

iBird et ([20121) noticed that the processes of stel¬ 
lar migration induced by accretion of a population of 
dark matter sub-halos are distinct from those in an iso¬ 
lated galactic disk. They found that the migration prob¬ 
ability of particles traces the radial mass distribution 
in the case of an isolated disk while in the perturbed 
disk this correlation is absent. Our simulations reveal 
that the migration probability of particles in a perturbed 
disk depends on the mass and orbital parameters of the 
merged satellite (Figure [12} . The migration probability 
of disk particles correlates with disk surface density be¬ 
yond i? ~ 7 kpc for the case of an isolated and some of 
the disturbed disks (see Figure [12] and top left panel of 
Figure[2]for a comparison with surface density). Such ra¬ 
dial profiles of the particle migration probability are sim¬ 
ilar to the one ass umed in the chemica l evolu tion model 
of our Galaxy of iSchonrich fc: BinnevI (I2009D . For the 
case of some other mass and orbital parameters of the 
merged satellite (e.g., when the satellite’s orbit lies in 
the disk plane) the migration probability of particles re- 
m ains about c onsta nt at i? ~ 7 to 20 kpc as it found 
by iBird et al.l (I2012D . but decreases with radius beyond 
R'^ 20 kpc. 

Figure [13] shows the angular momentum distribution 
evolution for models without satellite (AOO) and with 
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Figure 8. Radial oxygen abundance profiles at 4 Gyr for a disk galaxy without satellite (solid black line) and with satellite (dotted lines) 
after merging. The satellite masses are 10® Mq (mass ratio 1:70) (top left panel), 3 X 10® Mq (mass ratio 1:23) (top right panel), 6 X 10® 
Mq (mass ratio 1:12) (bottom left panel), and 9 X 10® Mq (mass ratio 1:8) (bottom right panel). 


satellites of different masses from 10® Mq to 9 x 10® Mq 
(A02, B02, C02, D02). Both axes are in kpc units due to 
the normalization of the angular momentum by the circu¬ 
lar velocity in the disks. All rows except for the last one 
represent the successive changes in angular momentum 
within the time interval of At = 0.5 Gyr. The last row 
represents the total cumulative changes of the angular 
momentum from 0 to 4 Gyr. Figure [14] shows the same 
as the previous figure but for the angular momentum dis¬ 
tribution evolution for models with a fixed satellite mass 
of 6 X 10® Mq and for different initial orbital inclinations 
z = 0°, 30°, 60°, and 90° (models GOl - G04). When con¬ 
sidering the angular momentum changes of disk particles 
after merging (bottom rows of Figured^ and Figure 114)) 
we find that the maximal changes in guiding radius Ai?g 
decrease outwards even for models that show a fairly 
constant value of migration probability within R 7 
to 20 kpc. Those migration patterns lead to the num¬ 
ber of particles migrating to the outer disk region being 
higher than the number of those migrating in the oppo¬ 
site direction. Regarding metallicity profile changes this 
means an increase of the mean metallicity in the outer 
disk region as is evident from our metallicity profile plots 
discussed above. The comparison of angular momentum 
changes in short time bins of At = 0.5 Gyr for isolated 


and perturbed disks (Figure |T3] and Figure [T4l) indicate 
that the accretion of a satellite may result in the sup¬ 
pression of further changes of angular momentum of the 
disk particles . This can be attributed to the suppression 
of non-axisymmetric patterns induced by the central bar 
after merging (see discussion below). 

Thus, the role of satellite accretion in inducing radial 
migration of the disk particles strongly depends on the 
orbit orientation and mass of the satellite. In general, the 
combined effects of spiral structure and bar patterns with 
satellite-induced migration mechanisms can produce very 
sophisticated migration patterns. 

4.4. Fourier analysis 

We calculate the discrete Fourier transform of the 
bulge-bdisk surface density. We divide the galactic disk 
in concentric rings and define the Fourier amplitude co¬ 
efficient from m = 0 to m = 4 in each ring. 

Figure [15] shows the radial dependence of the mean 
Fourier amplitudes Ai/Aq, A 2 /A 0 , A 3 /A 0 , and A4/A0 
from the bulge-l-disk surface density at a time bin near 
4 Gyr (averaged from 19 snapshots at time intervals from 
3.8 to 4.2 Gyr) for the model AOO (without satellite) 
and the models A02, B02, G02, and D02 (with satellites 
of different masses). The gray regions are estimates of 
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Figure 9. Radial profiles of the oxygen abundance’s standard deviation at 4 Gyr for a disk galaxy without satellite (solid black 

line) and with satellite (dotted lines) after merging. The satellite masses are 10® Mq (mass ratio 1:70) (top left panel), 3 X 10® Mq (mass 
ratio 1:23) (top right panel), 6 X 10® Mq (mass ratio 1:12) (bottom left panel), and 9 X 10® Mq (mass ratio 1:8) (bottom right panel). 
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Figure 10. Oxygen abundance distribution histogram at 4 Gyr 
for galactocentric radii of 10 kpc (top panel) and 20 kpc (bottom 
panel) for models that have the same satellite initial orbital incli¬ 
nation i = 30° and velocity V = 100 km s“^ but different masses. 
The vertical black line indicates the initial abundance at the galac¬ 
tocentric radii of 10 kpc (top panel) and 20 kpc (bottom panel). 


Figure 11. Oxygen abundance distribution histogram at 4 Gyr 
for galactocentric radii of 10 kpc (top panel) and 20 kpc (bottom 
panel) for models that have the same satellite initial orbital inclina¬ 
tion i = 30° and a mass of 6 X 10® Mq but different initial velocity. 
The vertical black line marks initial abundance at the galactocen¬ 
tric radii of 10 kpc (top panel) and 20 kpc (bottom panel). 
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Table 4 

Model runs with different satellite masses and initial positions and the resulting disk oxygen abundance gradients. 


Model 

Mass [x 10 ^Mq] 

R [kpc] 

i [deg] 

Alg(0/H),n 

Alg(0/H)inner 

Alg(0/H)outer 

AGO 

0 

- 

- 

-0.0453±0.0010 

-0.0387±0.0017 

-0.0554±0.0007 

AOl 

1 

40 

0 

-0.0375±0.0007 

-0.0303±0.0006 

-0.0466±0.0006 

A02 

1 

40 

30 

-0.0439±0.0010 

-0.0361±0.0012 

-0.0573±0.0007 

A03 

1 

40 

60 

-0.0449±0.0009 

-0.0387±0.0017 

-0.0535±0.0008 

A04 

1 

40 

90 

-0.0451±0.0009 

-0.0385±0.0018 

-0.0527±0.0007 

A05 

1 

40 

120 

-0.0454±0.0010 

-0.0386±0.0018 

-0.0538±0.0005 

A06 

1 

40 

150 

-0.0455±0.0010 

-0.0389±0.0017 

-0.0558±0.0008 

A07 

1 

40 

180 

-0.0434±0.0009 

-0.0371±0.0015 

-0.0533±0.0005 

BOl 

3 

40 

0 

-0.0315±0.0004 

-0.0274±0.0009 

-0.0353±0.0011 

B02 

3 

40 

30 

-0.0364±0.0009 

-0.0262±0.0008 

-0.0519±0.0008 

B03 

3 

40 

60 

-0.0400±0.0012 

-0.0257±0.0002 

-0.0603±0.0011 

B04 

3 

40 

90 

-0.0416±0.0011 

-0.0297±0.0006 

-0.0596±0.0005 

B05 

3 

40 

120 

-0.0419±0.0010 

-0.0349±0.0011 

-0.0592±0.0006 

B06 

3 

40 

150 

-0.0427±0.0010 

-0.0347±0.0012 

-0.0567±0.0007 

B07 

3 

40 

180 

-0.0365±0.0006 

-0.0314±0.0006 

-0.0466±0.0005 

COl 

6 

40 

0 

-0.0265±0.0002 

-0.0240±0.0003 

-0.0310±0.0005 

C02 

6 

40 

30 

-0.0326±0.0006 

-0.0265±0.0007 

-0.0426±0.0007 

C03 

6 

40 

60 

-0.0331±0.0009 

-0.0242±0.0004 

-0.0516±0.0016 

C04 

6 

40 

90 

-0.0385±0.0011 

-0.0270±0.0011 

-0.0576±0.0008 

C05 

6 

40 

120 

-0.0391±0.0009 

-0.0302±0.0007 

-0.0546±0.0013 

C06 

6 

40 

150 

-0.0411±0.0006 

-0.0363±0.0007 

-0.0490±0.0005 

C07 

6 

40 

180 

-0.0342±0.0004 

-0.0300±0.0002 

-0.0431±0.0005 

DOl 

9 

40 

0 

-0.0228±0.0004 

-0.0154±0.0006 

-0.0247±0.0010 

D02 

9 

40 

30 

-0.0302±0.0005 

-0.0243±0.0004 

-0.0362±0.0011 

DOS 

9 

40 

60 

-0.0316±0.0006 

-0.0256±0.0005 

-0.0429±0.0018 

D04 

9 

40 

90 

-0.0370±0.0011 

-0.0239±0.0005 

-0.0570±0.0009 

D05 

9 

40 

120 

-0.0380±0.0008 

-0.0302±0.0013 

-0.0505±0.0010 

D06 

9 

40 

150 

-0.0372±0.0004 

-0.0365±0.0004 

-0.0395±0.0015 

D07 

9 

40 

180 

-0.0325±0.0002 

-0.0303±0.0004 

-0.0361±0.0005 


Note. — First column: Model series number. Column 2: Satellite mass in units of solar masses. Column 3: Initial distance of the 
satellite from the center of the disk galaxy. Column 4; Initial orbital inclination of the satellite. Column 5: Oxygen abundance gradients 
from fitting radial abundance profile across the disk radius range of H = 0 to 25 kpc in dex kpc“^. Column 6; Oxygen abundance gradients 
from fitting radial abundance profiles in the inner disk (ij = 5 - 15 kpc) in dex kpc“^. Column 7: Oxygen abundance gradients from 
fitting radial abundance profiles in the outer disk (i? = 15 - 25 kpc) in dex kpc“^. The first row of this table shows the results for a host 
galaxy without a satellite. 


the significance limits of the Fourier amplitudes, i.e., the 
standard deviation of the values of the Fourier ampli¬ 
tudes obtained from 19 snapshots at the mentioned time 
intervals from 3.8 to 4.2 Gyr. Since the standard de¬ 
viations of the values of the Fourier amplitudes for the 
models AGO, A02, B02, C02, and D02 have comparable 
radial profiles, we plot only the maximal value of the 
standard deviation of the presented models. In all cases 
halo and satellite particles were not taken into account. 
All Fourier amplitudes are normalized by the main 
axisymmetric component Aq. High amplitudes of A 2 IAq 
in the central disk region clearly confirm the existence 
of the central bar. Our simulations show that the host 
galaxy bar strength (in the terms of their Fourier am¬ 
plitudes A 2 /A 0 ) becomes weaker for lower-mass ratios of 
the merged galaxies. 

5. CONCLUSIONS 

We investigate the dynamical interaction and merg¬ 
ing of a Milky Way-like galaxy (host galaxy) with low- 
mass satellites (mass ratio 1:70 - 1:8) during the last 
four Gyr. We use more than 7.4 million N-body parti¬ 
cles to study the dynamical evolution of such systems (2 
million particles for the disk, 0.4 million for the bulge, 
and 5 million for the halo of the host galaxy and between 
33000 to 300000 for a satellite). We examine the possible 


changes of the radial abundance gradient in the disk of 
the host galaxy caused by particle migration induced by 
the merger or interaction. The initial abundance distri¬ 
bution in the disk of the host galaxy is assumed to be an 
exponential (linear in log(0/H) scale). We computed a 
large set of mergers with different initial configurations, 
i.e., with different satellite masses, positions, and orbital 
velocities. We consider the evolution of a system over 
the period of 4 Gyr without star formation. 

We find that there are no significant metallicity 
changes at any radius due to particle migration in the 
case of the accretion of a low-mass satellite of 10 ® Mq 
( mass ratio 1:70) except for the special case of prograde 
satellite motion in the disk plane of the host galaxy. In 
general, the change in the abundance distribution in the 
disk of the host galaxy depends not only on the satellite 
mass but also on the initial inclination of the satellite or¬ 
bit. The largest changes take place if the satellite orbit 
lies in the disk plane of the host galaxy with prograde 
motion {i = 0°), i.e., the satellite moves through the 
disk in the direction of the host galaxy rotation. The 
smallest (if any) changes take place if the satellite has 
an initial retrograde orbital motion with an orbital in¬ 
clination i = 90° to 150°. Abundances at galactocentric 
distances larger than ^ 10 kpc can show a significant in¬ 
crease in the case of accretion of a satellite with a mass 
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Figure 12. The migration probability of disk and bulge particles. The lines correspond to the fraction of particles that change their 
galactocentric radius by more than 3 kpc after 4 Gyr of evolution. 


> 3 X 10® M 0 (mass ratio 1:23). 

The merger can result in a change of the slope of 
the abundance gradient. The radial abundance gradi¬ 
ent within the central region of the disk of the host is 
rather weakly affected by a merger with a satellite of a 
mass up to 9 X 10® Mq. The radial abundance gradient 
flattens in the range of galactocentric distances from 5 to 
15 kpc in the case of a merger with a satellite of a mass 

> 3 X 10® Mq. There is no significant change in the 
abundance gradient slope in the outer disk 15 kpc 
- 25 kpc) in any merger while the scatter in metallici- 
ties at a given radius increases significantly for most of 
the satellite initial masses/positions compared to the iso¬ 
lated galaxy case. Assum ing a mass-to-ligh t ratio in the 
B band of M/Lb = 1.4 (|Flvnn et al.ll200CT . the optical 
(i? 25 ) radius of our host galaxy is 15 kpc. This argues 
against attributing the break (flattening) of the abun¬ 
dance gradient near the optical radius observed in the 
extended disks of Milky Way-like galaxies exclusively to 
merger-induced stellar migration. 
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